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ABSTRACT 

We present three-integral dynamical models for the three Local Group dwarf elliptical 
galaxies: NGC147, NGC185, and NGC205. These models are fitted to the 2MASS J- 
band surface brightness distribution and the major-axis kinematics (mean streaming 
velocity and velocity dispersion) and, in the case of NGC205, also to the minor-axis 
kinematics. The kinematical information extends out to 2 in the case of NGC205 
and out to about 1 Rc in the case of NGC147 and NGC185. It is the first time 
models are constructed for the Local Group dEs that allow for the presence of dark 
matter at large radii and that are constrained by kinematics out to at least one half- 
light radius. The B-band mass-to-light ratios of all three galaxies are rather similar, 
(M/L)b « 3 - 4Mq/L0,b. Within the inner two half-fight radii, about 40 - 50% of 
the mass is in the form of dark matter, so dEs contain about as much dark matter as 
bright ellipticals. 

Based on their appreciable apparent flattening, we modeled NGC2G5 and NGC147 
as being viewed edge-on. For NGC185, having a much rounder appearance on the sky, 
we produced models for different inclinations. NGC205 and NGC147 have a relatively 
isotropic velocity dispersion tensor within the region where the internal dynamics are 
strongly constrained by the data. Our estimated inclination for NGC185 is i « 50° 
because in that case the model has an intrinsic flattening close to the peak of the 
intrinsic shape distribution of dEs and it, like the best fitting models for NGC147 and 
NGC205, is nearly isotropic. We also show that the dynamical properties of the bright 
nucleus of NGC205 are not unlike those of a massive globular cluster. 

Key v^fords: galaxies: individual: NGC147, NGC185, NGC205 - galaxies: dwarf 



1 INTRODUCTION 

From the days of the great comet-hunter Charles Messier up 
to well into the 20th century, NGC205 was the only known 
member of the class of objects that we now call dwarf ellip- 
ticals (dEs) (Messier 1801). dEs are faint (Ms > -18 mag) 
galaxies with diffuse, approximately exponentially declining 
surface brightness profiles and smooth elliptical isophotes 
(Ferguson & Binggeli 1994). NGC205 was resolved into stars 
and thus confirmed as a Local-Group member by Walter 
Baade, using the 100 inch telescope at Mount Wilson Ob- 
servatory (Baade 1944a). He correctly identified NGC205, at 
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a projected distance of only 37' from M31, as a diffuse stel- 
lar system with a stellar population very much like that of 
the Galactic globular clusters. NGC147 and NGC185 form 
a close pair, separated by only 58' on the sky (van den 
Bergh 1998) and located some 7° away from M31. Baade 
confirmed NGC147 and NGC185 as members of the Lo- 
cal Group by resolving them into individual stars (Baade 
1944b). No new Local Group dEs have been discovered since 
then while the sensus of the Local Group dwarf spheroidals 
(dSph), which are even fainter (absolute B-band magni- 
tude Mb > — 13 mag) and more diffuse than dEs, is most 
likely not yet complete. Since Baade could resolve these 
dwarf galaxies at about the same magnitude as M31, he 
surmised that they should be at roughly the same distance 
as M31. More precise distance measurements have since 
then clarified the geometry of the M31 satellite galaxy sys- 
tem. Using the TRGB distances measured by McConnachie 
(2005), -DNGC147 = 675 ± 27 kpc, Dngciss = 616 ± 26 kpc, 
I>NGC205 = 824 ± 27 kpc, Dm3i = 785 ± 25 kpc, which 
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Figure 1. 2MASS J-band image of NGC205. Overplotted onto 
tiiis image are the fitted elliptical isophotes at J-band surface 
brightness levels of 18, 19, 20, 21, 22, and 23 mag arcsec"^. Be- 
yond a major-axis distance of 300" the isophotes of NGC205 
start to twist, by almost 10° between the 21 to 23 mag arcsec"^ 
isophotes. The direction towards the center of M31 is indicated 
by an arrow. 



are in good agreement with those presented by Lee et al. 
(1993) and Lee (1996), and the angular distances between 
the galaxies, the following picture emerges. Using the dis- 
tance moduli and the sources of photometric errors given 
by McConnachie et al. (2005), we find that NGC147 and 
NGC185 arc separated by 63 ± 33 kpc and the couple is 
situated 160 kpc in front of M31. For all practical pur- 
poses, we can therefore treat them as isolated stellar sys- 
tems. NGC205, on the other hand, is situated 48 ± 30 kpc 
behind M31. 

In this paper, we present dynamical models for 
NGC147, NGC185, and NGC205. In order to construct a 
mass model for the luminous matter, we use 2MASS J-band 
images of these galaxies, downloaded from the NASA/IPAC 
Infrared Science Archive website. These dwarf galaxies have 
complex star-formation histories and show evidence for spa- 
tial variations of the stellar population, the most striking 
being the very blue nucleus of NGC205 (Butler &: Martmez- 
Delgado 2005; Han et al. 1997; Peletier 1993). Therefore, 
the density distribution obtained from near-infrared J-band 
images should trace the mass distribution of the luminous 
matter more closely than that obtained from optical images. 
We measured the surface-brightness profile, position angle, 
and cUipticity e = 1 — 6/a of all three galaxies as a function 
of the geometric mean of major and minor axis distance, 
denoted by a and b respectively. These were obtained using 
our own software. Basically, the code fits an ellipse through 
a set of positions where a given surface brightness level is 
reached. The shape of an isophote, relative to the best fitting 
ellipse with semi-major axis a and ellipticity e, is quantified 
by expanding the intensity variation along this ellipse in a 
fourth order Fourier series with coefficients S4, S3, Ca and 



C3: 

I{a, 9) = Io{a) [1 + C3{a) cos(36i)-F 

Ci{a) cos(46»)) + Saia) sin(36l)) + S^a) sin(46»)] . 



(1) 



Here, /o(o) is the average intensity of the isophote and the 
angle 6 is measured from the major axis. Using all photo- 
metric parameters, evaluated as a function of a, a smooth 
image of each galaxy is reconstructed. We used this smooth 
image to measure the total J-band magnitude and half-light 
radii. In each case, the total apparent magnitude we derive is 
within 0.05 mag of that listed in the 2MASS Large Galaxy 
Atlas (Jarrctt et al. 2003). As to the half-light radii, we 
find for NGC205 that i?o(NGC205) = 130" (0.52 kpc), for 
NGC147 that J?e(NGC147) = 122" (0.38 kpc), and NGC185 
that J?e(NGC185) = 90" (0.27 kpc). 

This smooth image is then deprojected in order to ob- 
tain the spatial density distribution of the stars. The spatial 
density of an axisymmetric stellar system cannot uniquely 
be reconstructed from the observed projected density if it is 
not viewed edge-on (Gerhard & Binney 1996). From the in- 
finity of equally plausible positive, axially symmetric spatial 
densities that project to the observed projected density, we 
pick the one that can be represented as a weighted sum of 
basis functions, p{zu,z) = '^^^Ck Pkii^, z), with pk{zu,z) of 
the form 



Pk{^,z) = exp 



(2) 



Here, is a scale radius, ruk is a shape parameter, mak- 
ing the density profile more (small mk) or less (large mk) 
centrally concentrated, q is the axis-ratio, controlling the 
fiattening of the basis function. These basis functions pro- 
duce elliptical isophotes when projected onto the sky. For 
a given inclination, the unknown coefficients in the expan- 
sion of the spatial density are determined by fitting the pro- 
jected basis functions to the observed photometry, evaluated 
on an elliptical, logarithmically spaced grid subject to the 
constraint that the spatial density is positive everywhere us- 
ing a Quadratic Programming algorithm (Dejonghe 1989). 
Note that the coefficients Cfc can be negative, as long as the 
weighted sum of haais functions is positive. This algorithm 
has a library at its disposal that contains hundreds of basis 
functions to choose from. Typically, the does not further 
improve after co-adding a few tens of components and the 
procedure can be stopped. This deprojection procedure en- 
sures that the spatial density is smooth, well-behaved, and 
positive. 

Already from the 2MASS image, it is obvious that 
NGC205 is noticeably influenced by the tidal forces exerted 
by M31 (see Fig. 1). Beyond a major-axis distance of « 300" 
the isophotes of NGC205 twist from a position angle of 167° 
to 160° at 400". At that radius, the surface brightness pro- 
file also shows a downward break, as can clearly be seen in 
Fig. 5 (see also Hodge (1973) and Choi et al. (2002)). In 
these outer regions, peculiar kinematics, such as the outer- 
most stars moving counter to the rotation of the main body 
(Geha et al. 2006), provide further evidence for a tidal in- 
teraction with M31. We will therefore check if our results 
for NGC205 are affected by it being so close to M31. The 
isophotes of NGC147 and NGC185 axe much more regular, 
with negligible isophote twisting, see Figs. 2 and 3. 
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Figure 2. 2MASS J-band image of NGC147. Ovcrplottcd onto 
tiiis image are the fitted elliptical isophotes at J-band surface 
brightness levels of 20, 21, 22, and 23 mag arcsec"^. 
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Figure 3. 2MASS J-band image of NGC185. Overplotted onto 
this image are the fitted elliptical isophotes at J-band surface 
brightness levels of 18, 19, 20, 21, 22, and 23 mag arcsec"^. 



The spectra were obtained with the 1.93-m telescope 
of the Observatoire de Haute-Provencc, equipped with the 
CARELEC long-sIit spectrograph. In combination with a 
EEV receptor with 2048 x 780 pixels of 15 /jm, the instru- 
mental velocity dispersion was ~ 25 km s"^ (R=5050) for 
a slit width of 1.5". More details and a log of the observa- 



tions can be found in Simien & Prugniel (2002). 10-minute 
exposures on blank fields of sky obtained during the same 
night as the science exposures were used for sky subtrac- 
tion. With a Fourier-Fitting technique (Franx et al. 1989), 
the spectrum of a G8111 template star (HD5395), convolved 
with a Gaussian line-of-sight velocity distribution, yielded a 
good fit to the galaxy spectrum in Fourier space, providing 
a simultaneous estimate of the mean streaming velocity Vp 
and the projected velocity dispersion (jp. 

From the deep major axis spectrum of NGC205, stel- 
lar kinematics out to about 300" (1.2 kpc or 2.3 half-light 
radii) could be extracted. The minor axis spectra yields use- 
ful kinematic information out to 120" (0.9 kpc). Thanks to 
the proximity of this galaxy and the unprecedented depth of 
the spectra, we can measure kinematics near the center of 
NGC205 with a spatial resolution of ~ 5 — 10 parsecs. The 
major axis spectrum of NGC185 allows extracting kinemat- 
ics out to ~ 100" (0.3 kpc or 1.1 R^). NGC147 has the lowest 
surface brightness of all Local Group dEs and is moreover 
riddled with foreground stars and was therefore the most 
difficult object to obtain decent kinematics for. The spectra 
of NGC147 secured by Simien & Prugniel (2002) turned out 
to be insufficiently deep. We, therefore, re-observed NGC147 
in Jarmary 2002 and January 2003, using the same instru- 
mental setup and the same data reduction procedure. With 
the addition of these new exposures (six one-hour exposures 
along the southern semi-major axis and four one-hour expo- 
sures along the northern semi- major axis), the major-axis 
kinematics of NGC147 now extend out to ~ 150" (0.5 kpc 
or 1.2 Re) along the northern semi-major axis and out to 
~ 100" (0.3 kpc or 0.8 Re) along the southern semi-major 
axis. 

The details of the modeling method are given in Section 
2. We present the results for each of the galax;ies in Sect. 3. 
Our conclusions are summarized in Sect. 4. 



2 THREE-INTEGRAL DYNAMICAL MODELS 

The internal dynamics of a steady-state oblate stellar sys- 
tem are described by a gravitational potential V'(w, z), that 
determines the stellar orbits, and the distribution function 
(DF) F{r, v) dfdv, which gives the number density of stars 
in phase space {{zu,z) arc cylindrical coordinates in each 
meridional plane). Loosely speaking, the DF distributes the 
stars over all possible orbits in a given potential. In the fol- 
lowing, we will work in prolate elliptical coordinates (A, j/) 
(Dejonghe & de Zeeuw 1988), defined by 



2 _ (A + (i)(;/ + q) 



(A + i,)('^ + 7) 



(3) 



a — 7 7 — a 

with —7 < v < —a < A. Surfaces of constant A are 
confocal prolate spheroids; surfaces defined by constant v 
are confocal two-sheeted hyperboloids, both with foci at 
z = ±A = ±V7 — «• We will approximate the gravitational 
potential by a Stackel potential, of the form 



ip{^,z) = 



(A + 7)G(A)-(z/ + 7)G(i/) 



(4) 



with G(A) the potential in the equatorial plane. Such po- 
tentials allow the existence of three integrals of motion: the 
binding energy E, the z-component of the angular moment 
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Lz, which we use in the more convenient form I2 = L'l/2, 
and 

/3 4(L^-L^.) + (7-«)(i..^-.^^^^^^), (5) 

an integral which is a generalisation of angular momentum 
conservation in spherical systems. Roughly, for a given E 
and I2, 13 determines how high above the equatorial plane an 
orbit will come. Models that make use of the third integral 
are called three-integral or 31 models, versus two-integral or 
21 models that go without it. 

A detailed account of the method we employed to con- 
struct the spheroidal coordinate system and Stackel poten- 
tial that give the best fit to a given a^cisymmetric potential 
can be found in Dejonghe & de Zeeuw (1988), Dejonghe et al. 
(1996), and in De Bruyne et al. (2001). In brief, we deproject 
the observed surface brightness distribution, derived from a 
2MASS J-band image, assuming the galaxy to be axisym- 
metric. Since NGC147 and NGC205 have a fairly flattened 
apparent shape, with an apparent axial ratio q^app = 0.5, 
it is quite unlikely that they are viewed far from edge-on. 
The intrinsic flattening distribution of dEs peaks around 
an intrinsic axial ratio gintr = 0.6 while galaxies with an 
intrinsic axial ratio smaller than gintr = 0.4 are virtually 
absent (Binggeli & Popescu 1995). Assuming both galaxies 
to have an oblate, axissymetric light distribution that is in- 
trinsically rounder than a gintr = 0-4 shape, the inclination 
is expected to be larger than i ^ 70°. Therefore, we can as- 
sume for simplicity that NGC147 and NGC205 are viewed 
edge-on without significantly affecting any of the results. 
NGC185, on the other hand, has a much rounder apparent 
shape, with an apparent axial ratio gapp = 0.8. If we as- 
sume that the galaxy has an intrinsic shape that is rounder 
than gintr = 0.4, the inclination is expected to be larger 
than i ~ 40° . In order to assess the influence of the intrinsic 
shape and the viewing angle of NGC185 on our results, we 
generated models with i = 90°, i = 60°, i = 55°, i = 50°, 
and i = 45° . However, with kinematics along the major axis 
only and lacking higher-order moments of the line-of-sight 
velocity distribution, the kinematical data do not allow to 
constrain the inclination with any degree of precision. 

The total mass density, including dark matter, is pa- 
rameterized as the spatial luminous mass density multiplied 
by a spatially varying mass-to-light ratio 

^iw,z) = A(l + B^zu^ + {z/qry (6) 

with the parameters A and B to be estimated from the data 
and q the axis ratio of the luminosity density distribution. 
Attempts to find intermediate-mass black holes in dEs, using 
either ground-based kinematics (Geha, Guhathakurta, van 
der Marel 2002) or high spatial resolution kinematical infor- 
mation obtained with HST (Valluri et al. 2005), have sofar 
failed. For this reason, and since we arc mostly interested in 
the total dark-matter content of these galaxies, we did not 
include a central black hole in our models. The gravitational 
potential is obtained by decomposing the total mass density 
in spherical harmonics. Finally, we determine the Stackel 
potential, i.e. the focal length A and the function G, that 
best fits this gravitational potential. The Stackel form of 
the potential is used only to calculate the third integral I3; 
elsewhere the potential derived from spherical harmonics is 



used. The potential is normalised such that i/'(r) = 1/r at 
large radius r, with r expressed in kiloparsccs. 

For a given potential, we wish to find the DF that best 
reproduces the kinematical information. As a first step, the 
DF is written as a weighted sum of basis functions, called 
"components" . Here, we will use components of the form 

F'''-^{E,h,h) = {E-Eo,ir{l2-Io.irit, 

if E> Eo,i, eih > Io,i > 
= 0, iiE< Eo,i or eil2 < Io,i (7) 

with <Ji, Ti, and pi integer numbers. iJo,, is a lower bound 
on the binding energy, which defines the outer boundary 
within which the component is non-zero. The parameter e, 
indicates whether stars are placed only on orbits with pos- 
itive angular momentum (e^ = -1-1) or only on orbits with 
negative angular momentum {ti = —1) or on both {ti = 0). 
We refer the reader to Appendix A for a discussion of the 
properties of these components. The DF then takes the form 
F(E,l2,h) = X]»CiJ'''''(S,/2,/3). The coefficients a are 
determined by minimizing the quantity 

^,^^j^ obs,-E,c.obs; y^ (8) 

with obs( an observed data point, ai the l-cr errorbar on that 
data point, and obs; the corresponding value calculated from 
the basis function 7?*'*=', subject to the constraint that the 
DF be positive everywhere in phase space. Note that the 
coefficients Cj can be negative, as long as the weighted sum 
of basis functions is positive. 

The projected velocity moments, /i„ {Q.) can be obtained 
by integrating suitably weighted combinations of the spatial 
velocity moments over a line of sight in a direction f2: 

Atn(n) = / v;F{E,l2,h)dvd^ 
Jn 

= / v;F''''{E,l2,l3)dvd4 (9) 

with Vp the line-of-sight velocity and ^ a coordinate along 
the line of sight. From the observed quantities, the pro- 
jected luminosity density Pp = lO^^'^'^®"''''^/^ '^ Lj,opc"^ 
with hj the observed surface brightness expressed in mag- 
nitude arcsec"^, the mean projected velocity Vp, and the 
projected velocity dispersion ap, one can easily construct ob- 
served values for the lowest-order velocity moments /xo = Pp, 
Ml = PpVp, and /i2 = pp{vp + Cp) along each observed line 
of sight and insert these values into eqn. (8), which is mini- 
mized using a Quadratic Programming algorithm (Dejonghe 
1989). The minimisation algorithm has a library at its dis- 
posal that contains hundreds of components to choose from. 
Since the algorithm itself determines which and how many 
components are used, this is a virtually non-parametric 
method of retrieving the distribution function. Typically, 
the does not further improve after co-adding a few tens of 
components and the procedure can be stopped. This finally 
yields the coefficients c, and consequently the distribution 
function that best reproduces the data for a given gravi- 
tational potential. If the basis functions F''"' are smooth 
and well-behaved functions of the integrals of motion, any 
combination of a finite number of basis functions will by con- 
struction also be smooth and well-behaved in phase space. 
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Figure 4. Contours of constant ^^r the models fitted to the 
NGC205 data as a function of the total mass within a 2 i?e radius, 
M{2Rg), and the scale parameter B of equation (6). The white 
dots mark the parameter values of the models we constructed. 
The contours indicate the 64%, 90%, 95%, and 99% confidence 
levels. At the 95% confidence level, the mass within a 2 i?e radius 
is M(2 _Rc) = 10.2^2:2 X 10^ Mq. For the scale parameter, we find 
B = l.llt f kpc-i, at the same confidence level. The best fitting 
model evidently requires a spatially varying mass-to-light ratio to 
reproduce the data. Models with a constant mass-to-light ratio, 
i.e. with B = kpc~^ are barely acceptable at the 95% confidence 
level. The vertical lines indicate the locus of models for which the 
mass-to-light ratio at 2 Rc equals 1, 1.5, 2, and 2.5 times the 
estimated stellar mass-to-light ratio ((M/L)j = 1.32, see section 
3.1.2). All acceptable models, even the ones with a zero B, have 
mass-to-light ratios well above that of a purely stellar system. 

It is well known that in a spherical geometry even com- 
plete kinematical knowledge does not suffice to uniquely con- 
strain both the gravitational potential and the distribution 
function. On the contrary : for each potential, there is a 
unique distribution function that reproduces the data (De- 
jonghe & Merritt 1992). That distribution function, how- 
ever, need not be positive everywhere in phase-space. The 
positivity constraint we are imposing on the distribution 
function therefore helps to constrain the gravitational po- 
tential. The degeneracy in the general axisymmetric case, 
however, can be expected to be a lot smaller than in the 
spherical case. For 21 models, the even part of the distribu- 
tion function can be derived directly from the augmented 
stellar density p{zu, ip) which, in the absence of dark matter, 
can in principle be constructed directly by deprojecting the 
observed surface brightness distribution. This is clearly not 
possible in the spherical case. Therefore, even if the kinemat- 
ical data are restricted to the three lowest order moments 
of the distribution function along the major (and sometimes 
minor) axis, wo can expect to be able to put meaningful 
constraints on the dark matter content and mass, a point 
already made in Dejonghe et al. (1996). 31 models have a 
larger freedom in distributing stars over all possible orbits 
than 21 models, widening somewhat the range of possible 
solutions that is consistent with the data. 

We fitted models to the data of each galaxy with various 
values for the parameters A and B in eqn. (6). The model 
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with the lowest x'^-value, which most closely reproduces the 
data, is retained as the best model. The uncertainties on 
the dynamical quantities presented in this paper are derived 
from the range of allowed models (all models within the 95% 
confidence level are deemed acceptable). 

3 RESULTS 
3.1 NGC205 

3.1.1 New and previous kinematical work 

The major and minor axis J-band photometry and kinemat- 
ics of NGC205 as measured by Simien & Prugniel (2002) are 
presented in Fig. 5. These data were derived from spectra 
with a 25 km s"'^ velocity resolution, using a Fourier-Fitting 
technique. Detailed simulations showed that template mis- 
match has little effect on the kinematical parameters (Kol- 
eva et al. 2006). There is a slight degeneracy between the 
velocity dispersion and the metallicity: a template of too 
high a metallicity will give too largo a dispersion. The tem- 
plate star we used has a metallicity [Fc/H]= —0.5 (Prugniel 
& Soubiran 2001). Our simulations show that for a galaxy 
with [Fe/H]= -1 and a template with [Fe/H]= -0.5, the 
dispersion can be overestimated by ~ 12% at most. How- 
ever, the fitting procedure includes an additive continuum 
which decreases the mismatch by an order of magnitude. 
Age mismatch, as is likely to be the case in the peculiar 
nucleus of NGC205, is not expected to bias the kinematics. 

NGC205 has a very steep central density cusp or nu- 
cleus (see panel a of Fig. 5). Within the inner 3 arc- 
seconds, the J-band surface brightness becomes almost 2 
mag arcsec"^ brighter. This sudden density jump is reflected 
in a corresponding drop of the velocity dispersion by about 
10 km s~^ towards the center; from Up = 30 ± 2 km s~^ 
at a radius of 10" down to a central velocity dispersion 
of cTp = 20 ± 1 km s~^ (see panels b and c of Fig. 5). 
The nucleus of NGC205 therefore appears to be a round 
and very dense but dynamically cold substructure. Out- 
side the nucleus, the velocity dispersion was found to rise 
outwardly from 30±2 km s^^ near the center up to about 
45±7 km s~^ at 300". The rotation velocity rises to a mgix- 
imum of 20±5 km s~^ at 200" and declines again beyond 
that radius. This was the first clear detection of rotation 
in this galaxy. There is evidence that the mean velocity re- 
verses sign at a radius of about 300". This feature agrees 
with other observations (Hodge 1973; Geha et al. 2006) that 
suggest that the tidal forces of M31 only appreciably affect 
the outer regions of NGC205, beyond 5'. Moreover, the ab- 
sence of minor-axis rotation (Simien & Prugniel 2002) and 
the presence of a very dense nucleus, which is likely to scat- 
ter stars off box-orbits, which axe the backbone of any slowly 
rotating triaxial mass-distribution, into loop-orbits (Evans 
& CoUett 1994; Merritt & Quinlan 1998) both argue against 
NGC205 having a strongly triaxial mass distribution, justi- 
fying our assumption of an oblate geometry for the construc- 
tion of dynamical models. 

Previous work on the kinematics of NGC205 has yielded 
contradictory results. Hold et al. (1990) present major axis 
velocity and velocity dispersion profiles out to about 1', mea- 
sured from spectra with a 83 km s~^ per pixel spectral reso- 
lution and using a G5 star spectrum as template. The mean 
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Figure 5. Fit to ttie kinematics of NGC205. Panel a: the J-band surface brightness, /ij, along major and minor axis; panel b: the 
major-axis velocity dispersion, ap; panel c: the minor-axis velocity dispersion; panel d: the mean streaming velocity along the major 
axis, Vp. The black dots are the data points, the grey curves correspond to the best fitting dynamical model. The top axis of each panel 
is labeled with radius expressed in kiloparsecs. The model was fitted to the surface brightness on a logarithmically spaced grid covering 
the whole face of the galaxy but for clarity we only plot the major and minor axis profiles in this figure. The model accurately reproduces 
the central surface brightness pealc (see insets in panel a), and the corresponding central drop in the velocity dispersion (see insets in 
panels b and c). The major-axis velocity dispersion appears to keep rising very slowly out to the last data point, at 5'. 




U5 (kpc) U5 (kpc) 



Figure 6. Contours for the radial velocity dispersion, (Ttc, the tangential velocity dispersion, cr^, the difference between the radial and 
the vertical velocity dispersion (7^ — (Tz, and the mean tangential velocity, v^, in the meridional plane of the best fit dynamical model 
for NGC205. All contours are labeled in km s~^. 
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(t!_(kpc) tij_(kpc) 

Figure 7. Panel (a): the distribution function of the best fit dynamical model for NGC205, plotted in turning-point space for loop orbits 
in the equatorial plane, characterized by Iz = 0. The pericenter distance is denoted by the apocenter distance by vj+. Radial orbits 
populate the ro_ = axis while circular orbits lie along the diagonal, -uj- has the same sign as the angular momentum Lz ■ The contour 
labels are expressed in units of log(L0 j pc~^ (50km s"^)^''). from —4.5 to 1.0 with a 0.5 interval. The orbital distribution within the 
equatorial plane is clearly very nearly isotropic, with only a very slight excess of stars with positive angular momentum L^, causing the 
rotation of NGC205. Panel b: to have a more detailed look at this slight excess of stars on positive angular momentum orbits, we plotted 
F{E, Lz > 0, 0) — F{E, Lz < 0, 0). The contour labels are expressed in units of log(LQ j pc^'* (50 km s^^)^"*). At small radii, there is 
almost exact counter-rotation, with stars on positive and negative angular momentum orbits almost cancelling any rotation. At larger 
radii, the number of stars on positive angular momentum orbits starts to dominate and at a radius of ~ 1 kpc, there are about twice as 
many stars on right-handed near-circular orbits as there are on left-handed ones. 



velocity dispersion is about 60-70 km s~^, with only a very 
slight decline towards the center, reaching about 50 km s~^ 
in the very center. No significant rotation was detected. 
These authors deduce a mass-to-light ratio of (M/L)b ~ 
7 Mq/Lq^b- Bender et al. (1991) on the other hand found 
the velocity dispersion to rise almost linearly with radius, 
from ~ 30 km s^^ near the center to ^ 70 km s^^ at one 
half-light radius. At the very center, the velocity dispersion 
was found to drop to 20 km s~^. Again, rotation turned out 
to be negligible. The kinematics were derived at a 30 km s~^ 
spectral resolution using a single GSlll stellar template spec- 
trum. These authors estimate the B-band mass-to-light ra- 
tio at (M/L)b = 8 ± 2M©/Lo,B. In Held et al. (1992), the 
mean value of the major-axis velocity dispersion of NGC205, 
(a) = 42 km/s, is used to arrive at (M/L)b ~ 3.5 Mq/Lq^b- 
Carter & Sadler (1990) on the other hand retrieved an al- 
most flat velocity dispersion out to 1', constant at about 
45 km s~^. Very near the center, the velocity dispersion 
was observed to drop to ~ 15 km s~^. These kinematics 
were measured using a composite stellar template template 
from a galaxy spectrum with a very high spectral resolution 
(5 km s~^). These authors estimate the mass-to-light ratio 
at (M/L)b ~ 9.4M0/L0,B with an uncertainty of a factor 
of 2. 

3.1.2 The dynamical mass and mass-to-light ratio 

From the range of models that is compatible with the data at 
the 95% confidence level, we estimate the mass within a 2 i?e 
radius sphere at M{2Re) = IO.2I2 2 x 10® M©, correspond- 
ing to a J-band mass-to-light ratio within the inner 2 Re of 
(M/L)j = 2.2+°'4 Mq/Lq^j or a B-band mass-to-light ratio 
of (M/L)b = 4..5tl l Mq/Lq,b, using a B-J= 2.397 color, 
corrected for Galactic reddening (de Vaucouleurs et al. 1991; 
Jarrett et al. 2003). The mean metallicity of NGC205 is quite 



well determined at [Fo/H] « —0.9, using various techniques 
(McConnachic ct al. 2005; Butler & Martmez-Delgado 2005; 
Richer et al. 1998). In the central regions of our dynamical 
models, mass follows light. If the central mass-to-light ratio 
of the best fit model, (M/L)j = 1.32 Mq/Lq^j, is indicative 
of that of the stellar population and if wc adopt a mean 
metallicity [Fe/H] ~ —0.9, we estimate the mean age of 
the stellar population at ~ 9 Gyr (Worthey 1994), which 
is not unreasonable. The model with the lowest mass that 
is still consistent with the data at the 95% confidence level 
has a constant mass-to-light ratio (M/L)j — 1.8 Mq / Lq^j , 
which is significantly higher than that estimated for the stel- 
lar population above. Hence, we can state that NGC205 is 
not a pure stellar system, completely devoid of dark mat- 
ter. Within a 2 Re sphere, NGC205 consists of about 60% 
luminous matter, by mass, and of 40% dark matter. 

3.1.3 The internal dynamics 

The best fitting 31 model for NGG205 very accurately re- 
produces the photometric and kinematical data, as can be 
judged from Fig. 5, where the kinematics of the best fitting 
model are overplotted onto the data. Since the dynamics 
of this galaxy are constrained by kinematical information 
along both the major and minor axes, it is instructive to 
more closely examine some internal dynamical quantities. 
In Fig. 6, the spatial kinematics arc presented in the merid- 
ional plane: the radial velocity dispersion , the tangential 
velocity dispersion cr^, the difference between the radial and 
vertical velocity dispersions cr^ — Oz (in the case of a 21 
model: a-cj — Oz = 0) , and the mean tangential velocity . 

Geometrically, it is clear that the minor-axis velocity 
dispersion puts strong constraints on the behaviour of the 
radial velocity dispersion along the symmetry axis while the 
major-axis kinematics constrain both the radial and the tan- 
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Figure 8. The internal dynamics of the best fit model for 
NGC205, plotted as a function of radius within the equatorial 
plane. Top panel; the logarithm of the luminosity density of the 
stars, p; middle panel: the rotation velocity, u^; bottom panel: the 
radial, tangential, and vertical components of the velocity disper- 
sion tensor. 



gential velocity dispersion within the equatorial plane. The 
vertical velocity dispersion is only very weakly constrained 
by the kinematical data; it will mostly be determined by 
the photometry, i.e. by the flattening of the stellar distribu- 
tion. As is clear from Figs. 5, 6, and 8, the behaviour of the 
major-a:xis projected velocity is reflected in that of the tan- 
gential component a^, of the velocity dispersion tensor: the 
maximum around 0.5 kpc (125") and the ourward rise be- 
yond 0.8 kpc (200") can be clearly identified. The run of 
the radial velocity dispersion along the 2;-axis closely traces 
the minor-axis projected velocity dispersion, as expected; 
cf. the maximum around 0.2 kpc (50"). In a 31 model, the 
radial and vertical velocity dispersions are decoupled. Still, 
the radial versus vertical anistropy is very small : the ver- 
tical velocity dispersion is only marginally smaller than the 
radial dispersion within the inner 0.8 kpc (200"). If we use 
the anisotropy parameter /J^.z = 1 — (f^/cfro)^ to quantify 



this anisotropy, then we find that —0.6 < /^^.z < 0.3 every- 
where in the meridional plane. The radial versus tangential 
anisotropy is somewhat larger. If we introduce the quantity 
/3ro,<#. = 1 — (o-^/cr^)^, then we find that —1 < /Jiu,^ < 0.1 
within the region that is well constrained by the kinematical 
data. Beyond a radius of 0.8 kpc, the tangential anisotropy 
increases to keep up with the observed rise of the projected 
velocity dispersion along the major axis. 

The mean tangential velocity rises only very slowly 
as a function of radius, as does the projected mean velocity 
Vp along the major cixis. The velocity reaches a maximum 
at 0.8 kpc (200") and declines beyond that radius (see Fig. 
8). Extrapolating somewhat beyond the extent of our data, 
the velocity is expected to reverse sign at about 300", in 
agreement with Geha (2006). From Fig. 7, it is clear that 
there is only a slight excess of stars with positive angular 
momentum causing the observed rotation. Only beyond 
~ 0.7 kpc docs the number of stars with > clearly 
dominate that of Lz < stars. At a radius of ~ 1 kpc, there 
arc about twice as many stars on right-handed near-circular 
orbits as there are on left-handed ones. 

NGC205 has a very distinct nucleus, which shows up in 
the dynamical model as a spherically symmetric, dynami- 
cally very cold substructure (see Figs. 6 and 8). By singling 
out those basis functions for the distribution function that 
make up the nucleus, i.e. the basis functions with a very 
small outer boundary, we can study the nucleus in more de- 
tail. The nucleus can be modeled as an isotropic star cluster 
with a one-dimensional velocity dispersion cr « 7 km s~^ 
and a truncation radius of about 50 pc. By integrating the 
spherically symmetric mass density of the nucleus, we es- 
timate its total mass at 1.4 x 10® Mq. This is in excellent 
agreement with a previous estimate based on a King-model 
fit to the surface brightness profile of the nucleus and its pro- 
jected velocity dispersion (Carter & Sadler 1990; Jones et al. 
1996). It should be noted, however, that the mass-to-light 
ratio of the nucleus is probably different from that of the 
main body of NGC205 (VaUuri et al. 2005). The functional 
form for the mass-to-light ratio adopted by us, eq. (6), can- 
not reproduce a central M/L- variation. However, this is not 
our intention: we are interested in quantifying the amount 
of dark matter at large radii. Since the nucleus contributes 
only of the order of 0.1 % of the total mass of the galaxy, 
our total mciss estimates are robust. 

Thus, the nucleus of NGC205 is structurally and dy- 
namically not unlike a massive globular cluster. The orbital 
decay time scale for a globular cluster in NGC205 by dy- 
namical friction is of the order 10* — 10^ years, significantly 
less than a Hubble time (Lotz et al. 2001; Jones et al. 1996; 
Tremaine, Ostriker, Spitzer 1975), which makes it plausible 
that one or more star clusters have spiraled inwards and 
settled at the bottom of the gravitational well. Though de- 
tailed N-body simulations of this process are still lacking, 
the observed properties of the NGC205 nucleus are not in- 
consistent with this scenario. 



3.1.4 The tidal influence of MSI 

The TRGB distance estimates from McConnachie (2005) 
place NGC205 48 ± 33 kpc behind M31. The la error on 
the TRGB estimates for the distances of NGC205 and M31 
are large enough that it is quite plausible that NGC205 and 
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Figure 9. The gravitational potential well V oi NGC205 in the 
equatorial plane as a function of radius zu. The top curve is 
the gravitational potential of the best fit dynamical model for 

NGC205 in the abscence of M31. From top to bottom, the other 
curves indicate the "effective" potential in the jjresence of the ax- 
isymmetrized tidal field of M31 (eq. (10)) for relative distances of 
40 (which essentially overlaps with the unperturbed potential of 
NGC205), 20, 10, and 5 kpc between NGC205 and M31. 

M31 arc at the same distance from the Milky Way and hence 
at a relative distance of a mere 10 kpc. 

In order to get a rough estimate of how the presence 
of M31 affects our results, we approximate the gravitational 
potential of M31 as that of a spherical mass distribution at 
a distance R from NGC205. For computational ease, we will 
further assume that the center of mass of M31 lies within 
the equatorial plane of NGC205 and that distances within 
the stellar body of NGC205 arc much smaller than R. We 
then axisymmetrize this potential by averaging it over the 
azimuthal angle 0. We finally obtain the expression 

for the axisymmetric potential corresponding to the tidal 

field of M31, up to fourth order in z/R and vj/R, with 
M{R) the mass of M31 enclosed within a sphere of radius 
R (Klypin, Zhao, Somcrvillc 2002). 

We plot the gravitational potential of the best fit dy- 
namical model for NGC205 as a function of radius in the 
equatorial plane. The top curve is the potential of the model 
in the absence of M31. From top to bottom, the other curves 
indicate the potential in the presence of the axisyrnmetrized 
tidal field of M31 for relative distances of 40, 20, 10 and 5 kpc 
between NGC205 and M31. Even for the minimum current 
relative distance allowed by the observations, i.e. 10 kpc, the 
tidal influence of the giant spiral galajcy is negligible within 
the region for which we have kinematical data, that is within 
the inner 1.2 kpc. In that case, stars beyond 2.8 kpc (700") 



are unbound. If NGC205 has passed by M31 with a pericen- 
tric distance as small as 5 kpc, stars are unbound beyond a 
radius of ~ 1.5 kpc. These values can be compared with the 
analytical estimate 

with Tperi the orbital radius of a satellite galaxy with mass 
Mgat around a gala:xy with mass Mgai, and a; = 9 or 1, for 
stars inside the satellite galaxy on prograde or retrograde 
circular orbits, respectively (Read ot al. 2006). For a peri- 
centric distance of 10 kpc, this yields rtidai ~ 1-0 to 2.0 kpc 
(depending on x); for a pericentric distance of 5 kpc, this 
yields nidai = 0.6 to 1.3 kpc. We substituted the mass of 
M31 enclosed within a radius rperi for Mgai. These are also 
the radii whore in Fig. 9 the "effective" potential starts to 
deviate significantly from that of NGC205 alone. This makes 
us confident that we arc providing an adequate description 
of the tidal enfiuence of M31 on NGC205. 

Small pericentric distances, i.e. smaller than ~ 5 — 
10 kpc, which may have occurred in the past, lead to a large 
effect on the force field within the inner 1.2 kpc but also 
give rise to a tidal radius well within the region constrained 
by the photometric and kinematical data, in which case the 
construnction of equilibrium models itself becomes mean- 
ingless. In short, our results are robust if NGC205 has kept 
a distance of at least ~ 10 kpc from M31. 

3.2 NGC147 

3.2.1 New and previous kinematical work 

The major and minor axis J-band photometry and the 
major-axis kinematics of NGC147 as measured by Simien 
& Prugniel (2002), including the additional exposures ob- 
tained in January 2002 and 2003, are presented in Fig. 10. 
The kinematics of NGC147 extend out to ~ 150" (0.5 kpc 
or 0.8 Re). The projected velocity dispersion is roughly con- 
stant, at about ap ~ 25 km s^^. The rotation curve is consis- 
tent with NGC147 having zero rotation. The velocity disper- 
sion profile is in good agreement with the one presented by 
Bender ct al. (1991). However, these authors found signifi- 
cant rotation in NGC147, with an amplitude of ~ 10 km 
at a radius of 100". They estimated the B-band mass-to-light 
ratio of NGC147 at {M/L)b =7 ±3 Mo/L©,b. 

3.2.2 The dynamical mass and mass-to-light ratio 

The total mass within a 2 _Ro radius sphere is esti- 
mated at M(2 7?e) = S.Otl'l ^ 10* M0, corresponding to 
a J-band mass-to-light ratio within the inner 2 Re of 
(M/L)j = 3.4l2 g Af0/L0,j or a B-band mass-to-light ra- 
tio of (M/L)b = 4.01^4 Af0/L0,B. The mean metallic- 
ity of the stellar population is [Fo/H] ~ —1.0 (Han et al. 
1997; McConnachie et al. 2005). There is no evidence that 
NGC147 contains a population of intermediate-age stars 
(Renzini 1998); Mould, Kristian, and Da Costa (1983) es- 
timate that 90% of the stellar population is older than 
12 Gyr. Using 12 Gyr as a rough estimate for the mean 
age of the stellar population, the stellar mass-to-light rar- 
tio is (M/L)j ~ 1.6Mq/L0,j or (M/L)b ~ 2.6Mq/Lq,b 
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Figure 11. The radial, tangential, and vertical components of the 
velocity dispersion tensor of the best fit model of NGC147. All 
quantities are plotted as a function of radius within the equato- 
rial plane. The spatial stellar density profile within the equatorial 

plane is consistent with being exponentially declining from a cen- 
tral density of 0.32 Lj_q pc~'^ with an exponential scale-length of 
0.31 kpc. 
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Figure 10. Fit to the kinematics of NGC147. Panel a: the J- 
band surfa<;e brightness, ixj, along major and minor axis; panel 
b: the mean streaming velocity along the major axis, Vp; panel 
c: the major-axis velocity dispersion, ap. The black dots are the 

data points, the grey curves correspond to the best fitting dy- 
namical model. The model was fitted to the surface brightness 
on a grid covering the whole face of the galaxy but for clarity 
we only plot the major and minor axis profiles in this figure. The 
major-axis velocity dispersion remains approximately constant at 
<Tp ~ 25 km out to the last data point, at 2.3'. 



3.3 NGC185 

3.3.1 New and previous kinematical work 

The major and minor axis J-band photometry and major- 
axis kinematics of NGC185 as measured by Simien & Prug- 
niel (2002) are presented in Fig. 12. The kinematics of 
NGC185 extend out to ~ 100" (0.3 kpc or 1.1 7?e). The 
major-cixis velocity dispersion is roughly constant at about 
30 km s~^, in good agreement with Bender et al. (1991). 
There is a hint of rotation, at variance with Bender ct al. 
(1991), but, as in the case of NGC147, sky subtraction at 
these low surface-brightness levels is very difficult and this 
might affect the measurements. Bender et al. (1991) estimate 
the B-band mass-to-light ratio of NGC185 at {M/L)b = 
5 ± 2Mq/L0,b. Held ct al. (1992) found an essentially flat 
velocity dispersion profile, with a ^ 30 km/s. Interest- 
ingly, their major-axis velocity dispersion profile shows a 
central peat, with the dispersion rising to a central value 
of 50 km/s. These authors used a spherically symmetric, 
isotropic model with an exponentially declining density pro- 
file to estimate the central mass-to-light ratio of NGC185 at 
{M/L)b^3Mq/Lq,b. 



(Worthey 1994). Comparing this with the dynamical mass- 
to-light ratio of the best fit dynamical model, we find that 
NGC147 consists of about 50% luminous matter, by mass, 
and of 50% dark matter. 



3.2.3 The internal dynamics 

The stellar density and the radial, tangential, and vertical 
components of the velocity dispersion tensor of the best fit 
model for NGC147 arc presented in Fig. 11. The velocity dis- 
persion tensor is almost isotropic, with ct ~ 20 — 25 km s^^, 
up to a radius of 0.3 kpc, with the vertical velocity disper- 
sion being somewhat lower to help flatten the galaxy. Out- 
side 0.3 kpc, the tangential velocity dispersion rises up to 



3.3.2 The dynamical mass and mass-to-light ratio 

The best fitting model, with i = 50° and an intrinsic axial 
ratio gintr = 0.62, yields our best estimate for the mass of 
NGC185. The error-bar on this quantity is estimated taking 
into account all inclination angles. We estimate the mass 
within a 2 _Rc radius sphere at M{2Rc) = 2.6to[l x 10® M©, 
corresponding to a J-band mass-to-light ratio within the in- 
ner 2 Re of (M/L)j = 2.2toi Mq/Lq.j or a B-band mass- 
to-light ratio of (M/L)b = 3.0to:? Mo/Lq.b. This result 
is in good agreement with the one obtained by Held et al. 
(1992). The mean mctallicity of the stellar population is es- 
timated at [Fe/H] = -1.43 ± 0.15 by Martmez-Delgado & 
D. & Aparicio (1998). Lee et al. (1993) find NGC185 to 
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Figure 12. Fit to the streaming velocity along the major axis, Vp, 
the major and minor axis velocity dispersion, Up, and the J-band 
surface brightness, fxj, of NGC185 with a 31 dynamical model 
with inclination i = 50°. The black dots are the data points, the 
grey curves correspond to the model. The model was fitted to the 
surface brightness on a grid covering the whole face of the galaxy 
but for clarity we only plot the major and minor axis profiles in 
this figure. 



contain a complex stellar population with a mean metallic- 

ity of [Fc/H] = -1.23 ± 0.16. McConriachic ct al. (2005) 
derive a value [Fe/H] = —1.2. Using 10 Gyr as a rough esti- 
mate for the mean age of the stellar population (Martmez- 
Delgado, Aparicio, Gallart 1999), the stellar mass-to-light 
ratio is (M/L)j ~ l.SMo/i©,,! or (M / L)b ~ 2.1Mq/Lq,b 
(Worthcy 1994). In other words, NGC185 consists of about 
60% luminous matter, by mass, and of 40% dark matter. 



5.5.5 The internal dynamics 

The stellar density and the radial, tangential, and vertical 
components of the velocity dispersion tensor of the best fit 
models for NGC185 are presented in Fig. 13. The velocity 
dispersion tensor of the edge-on model {i = 90°, see left 
column of Fig. 13) has a remarkably large vertical velocity 



dispersion, peaking up to Uz ~ 40 km s^^, in the inner parts 
of the galaxy. This feature is not just an oddity of this one 
dynamical model but is common to all edge-on models we 
fitted to NGC185. Since within the context of axisymmet- 
ric models, viewed edge-on, the vertical velocity dispersion 
is determined almost solely by the photometry, we have to 
conclude that if NGC185 is viewed edge-on, it can only main- 
tain its round shape by having the vertical velocity disper- 
sion much larger than the radial velocity dispersion, which 
is set by the observed major-axis velocity dispersion. The 
model with i = 50° and an intrinsic axial ratio gintr = 0.62, 
is the most nearly isotropic of all models we constructed 
for NGC185, sec the middle column of Fig. 13. The rota- 
tion velocity of the stars flattens at about 22 km s^^ in this 
model, versus a rotation velocity of 8 km s~^ in the edge-on 
model. Since (i) the intrinsic axial ratio of this model is still 
fairly close to the peak of the intrinsic shape distribution of 
dEs, and {ii), like the best fitting models for NGC147 and 
NGC205, it is very nearly isotropic, this is our preferred best 
model for NGC185. The intrinsically very flattened model, 
with 9intr = 0.5 and i = 45° , see the right column of Fig. 12, 
has a much larger rotation velocity, its peak projected ro- 
tation velocity is somewhat higher than that of the i — 90° 
model but still safely within the error-bars. The vertical ve- 
locity dispersion CTz is now much smaller than in the less 
inclined models. Although this model agrees with the obser- 
vations to within the (quite substantial) error-bars, it fails to 
reproduce the outward rise of the major-axis velocity disper- 
sion. Hence, with better data it should be possible to better 
constrain the inclination of NGC185 on purely dynamical 
grounds. 



4 CONCLUSIONS 

In this paper, we present dynamical models for NGC147, 
NGC185, and NGC205. NGC205 and NGC147 have already 
a flattened apparent shape so we modeled them assuming 
them to be viewed edge-on [i = 90°). NGC185, on the other 
hand, has a much rounder projected shape than the two 
other Local Group dEs, so we produced models for it as- 
suming inclination angles i = 90° , 60° , 55° , 50° , and 45° . 

Ftom the range of models that is compatible with 
the data at the 95% confidence level, we estimate 
the mass of NGC205 within a. 2 Re radius sphere 



at M{2Re 



IO.2I2 Q X 10* M0, corresponding to a 



J-band mass-to-light ratio within the inner 2 Re of 
(M/L)j = 2.2+°-l Mq/Lo.j or a B-band mass-to-light ratio 
of (M/L)b = 4.5tl-l Mq/Lq,b- The best fitting 31 model for 
NGC205 very accurately reproduces the photometric and 
kinematical data, such as the very bright nucleus and the 
central drop in the velocity dispersion. Overall, the velocity 
dispersion tensor is quite isotropic. Only beyond a radius of 
0.8 kpc docs the tangential anisotropy increase to keep up 
with the observed rise of the projected velocity dispersion 
along the major axis. The very bright nucleus shows up in 
the model as spherically symmetric, dynamically very cold 
substructure, which can be modeled eis an isotropic star clus- 
ter with a total mass of 1.36 x 10^ M©, a one-dimensional 
velocity dispersion u « 7 km s~^ and a truncation radius of 
about 50 pc. Thus, the nucleus of NGC205 is structurally 
and dynamically not unlike a massive globular cluster. 
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We estimate the total mass of NGC147 within a 2 _Rc 
radius sphere at M{2Re) = 2,.Qt\^ x 10* Mq, correspond- 
ing to a J-band mass-to-light ratio within the inner 2 Re 
of (M/L)j = 3.4110 M0/L0, J or a B-band mass-to-light 
ratio of (M/L)b = 4.0l2;4 Mc5/Lc?,b. As in the case of 
NGC205, the velocity dispersion tensor of the best fit model 
for NGC147 is remarkably isotropic. 

We estimate the total mass of NGC185 within a 2 7?c 
radius sphere at M(2 R^) = 2.61°;® x 10* Mq, corresponding 
to a J-band mass-to-light ratio within the inner 2 _Rc of 
{M/L)j = 2.2to 5 Mq/Lq.j or a B-band mass-to-light ratio 
of {M/L)b = 3^OlJ:?M0/L0,B. The model with i = 50° 
has an intrinsic flattening that is still close to the peak of 
the intrinsic shape distribution of dEs and it, like the best 
fitting models for NGC147 and NGC205, is nearly isotropic. 
Therefore, this is our preferred model for NGC185. 

Summarizing, our dynamical models for the three Lo- 
cal Group dEs show them to have nearly isotropic velocity 
dispersion tensors. They have very similar mass-to-light ra- 
tios, of the order of (M/L)b ~ A^Mq/Lq^^. While this is 
still larger than the expected mass-to-light ratio of the stel- 
lar populations of these galaxies, it is about a factor of two 
smaller than the mass-to-light ratios derived by many pre- 
vious studies. This rather large difference is due in part to 
differences in the data analysis. Depending on the spectral 
resolution, the signal-to-noise ratio of the spectra, and the 
method of extracting kinematical data from the spectra, the 
velocity dispersion estimates for the same galaxy can differ 
by a factor of 1.5 — 2 between different authors, with lower 
resolution observations tending to yield higher values for 
the velocity dispersion. This, exacerbated by the difficulty 
of measuring core-radii and surface brightnesses of such dif- 
fuse objects, required for use in the standard King formula 
(Richstone & Tremaine 1986), seems to account for at least 
part of the large spread on the M/L estimates. Therefore, 
these dwarf galaxies contain much less dark matter than was 
believed to be the case up to now. Still, within the inner two 
half-light radii, about 40 — 50% of the mass is in the form of 
dark matter, so dEs contain as much dark matter as bright 
ellipticals. Indeed, Gerhard et al. (2001) produced dynami- 
cal models for a sample of 21 luminous ellipticals and found 
that the dark matter contributes 10—40% of the mass within 
1 Re and equal interior mass of dark and luminous matter 
occurs at 2 — 4i?e. 
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APPENDIX A: THE THREE-INTEGRAL 
COMPONENTS 

For a given potential, we wish to find the DF that best 
reproduces the kinematical information. As a first step, the 
DF is written as a weighted sum of basis functions, called 
"components" . Here, we use components of the form 

ii E> Eo,i, eih > /o,« > 
= 0, if£:<Bo,iOrei72</o,i(Al) 

with Oi, Ti, and pi integer numbers. Eo,i is a lower bound 
on the binding energy, which defines the outer boundary 
within which the component is non-zero. The condition 
£1/2 > Io,i > 0, with ei € [—1,-1-1] defines the "handed- 
ness" of the component. If — +1, only orbits with posi- 
tive angular momentum I2 > lo arc populated and we call 
the component _F''+ "right-handed". If f,: = —1, only or- 
bits with negative angular momentum I2 < —lo are pop- 
ulated and the component F^'~ is "left-handed". Compo- 
nents that populate progradc and retrograde orbits equally, 
labeled with ti = 0, arc constructed as F''"{E, 12,13) = 
F'-+{E,l2,h) + F''-{E,l2,h)- Components with /o 7^ 
and Eo =^ populate orbits within a torus-shaped region of 
space. 

The spatial velocity moments, /"2;^„ 2n('^' 
component F^'^ {E, I2, 13) are given by: 



F^'^{E,l2j3)v\V^vT dvxV4,v^ 

2! + n+(m + l)/2 /■V' 



I- 



(A2) 



{E-Eo,iY'dE 



X / I^' 



j(m-l)/2 



{h-h,iY*dl2 



{1+ - l3)'-''\l3 - I^r-"^dl3, (A3) 



with 



(A + ,)(G(A)-i.-^) 



(A4) 

I3 = (^ + 7)(gH-S-^). (A5) 

This integration can be performed analytically by repeatedly 

using the lemma 

b 



I 



x\b - xy{x - aydx = B{k -h 1, j + 1) 



X {h-ay+''+''V2Fi[-i,3 + l-3 + k + 2-l 



,(A6) 



with B the Euler beta function and 2F1 the hypergeometric 
function. Thus, one obtains the following expression for the 
spatial velocity moments: 



M2itm,2n(w.2) 



iTi+Ti+n+l+2 



X [ip- £o,i)''"-'^/' [(A + 7)(G(A) - EoM 



(-pi)rG + l/2)r 



r!(o-j + Ti + n 



+ / + 3), 



(A-t^)(V>-go,, - Jo,»/ro^) 
(A-h7)(G(A)-So,i) 



(r — pi)s{ai + 1)b f tp — Eo,i — Ioa/ 



^ s\i(Ji+Ti + n + l + r + 



s=0 



3). V ^(A) - Eo,i 



X E 



p=0 



p\ \{a-'y){G{X)-Eo,i) 



X 2 



Fi (-P - ^^^^Y^, o-i + Ti + n + l + r + s + 2; 



(Ji + n + n + I + r + s + 3;1 



.(A8) 



wmj-Eo,i) 

with the forefactor J^lrn,n given by 

-pi _ 

l,m,n 

2„+z+(^+i)/2 r(a. + i)r(r. + i)r(n + |)ra + i) 

r(o-i + Ti + n + l + 3) 

The velocity moments of the left-liaudcd component 
F^'~{E, 12,13), fJ'2i~m2n^ linked to those of its right- 
handed analog ^''+(7?, 72, /s) via the relation H2i~m.2n = 
(— l)'"/i2i^^ 2n- 'I'h^ spatial velocity moments with respect 
to the cylindrical velocity components v^^, V4,, and Vz can 
be constructed from the moments /ijf j„ 2n using the trans- 
formation 



(A7) 



sign(2;) Vz 



cos — sin \ / v\ 



sin cos & I \ Vi, 



(A9) 



with 



cos 6 = i I -, -i^, sm6 = ^ / 7 ^ ^.(AlO) 



(a-7)(A- z/)' 



(7-a)(A-i^)' 



This paper has been produced using the Royal Astronomical 
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